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Abstract 

This paper considers a spectrum sharing based cognitive radio (CR) communication system, which consists 
of a secondary user (SU) having muhiple transmit antennas and a single receive antenna and a primary user 
(PU) having a single receive antenna. The channel state information (CSI) on the link of the SU is assumed 
to be perfectly known at the SU transmitter (SU-Tx). However, due to loose cooperation between the SU 
and the PU, only partial CSI of the link between the SU-Tx and the PU is available at the SU-Tx. With the 
partial CSI and a prescribed transmit power constraint, our design objective is to determine the transmit signal 
covariance matrix that maximizes the rate of the SU while keeping the interference power to the PU below 
a threshold for all the possible channel realization within an uncertainty set. This problem, termed the robust 
cognitive beamforming problem, can be naturally formulated as a semi-infinite programming (SIP) problem with 
infinitely many constraints. This problem is first transformed into the second order cone programming (SOCP) 
problem and then solved via a standard interior point algorithm. Then, an analytical solution with much reduced 
complexity is developed from a geometric perspective. It is shown that both algorithms obtain the same optimal 
solution. Simulation examples are presented to validate the effectiveness of the proposed algorithms. 
Keywords: Cognitive radio, interference constraint, multiple-input single-output (MISO), partial channel state 
information, power allocation, rate maximization. 

I. INTRODUCTION 

One of the fundamental challenges faced by the wireless communication industry is how to meet 
rapidly growing demands for wireless services and applications with limited radio spectrum. Cognitive 
radio (CR) technology has been proposed as a promising solution to tackle such a challenge [l]-[8]. 
In a spectrum sharing based CR network, the secondary users (SUs) are allowed to coexist with the 
primary user (PU), subject to the constraint, namely the interference constraint, that the interference 
power from the SU to the PU is less than an acceptable value. Evidently, the purpose of the imposed 
interference constraint is to ensure that the quality of service (QoS) of the PU is not degraded due to the 
SUs. To be aware of whether the interference constraint is satisfied, the SUs needs obtain knowledge 
of the radio environment cognitively. 

In this paper, we consider a spectrum sharing based CR communication scenario, in which the SU 
uses a multiple-input single-output (MISO) channel and the primary user (PU) has one receive antenna. 
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We assume that the channel state information (CSI) about the SU link is perfectly known at the SU 
transmitter (SU-Tx). However, owing to loose cooperation between the SU and the PU, only the mean 
and covariance of the channel between the SU-Tx and the PU is available at the SU-Tx. With this 
CSI, our design objective is, for a given transmit power constraint, to determine the transmit signal 
covariance matrix that maximizes the rate of the SU while keeping the interference power to the PU 
below a threshold for all the possible channel realizations within an uncertainty set. We term this design 
problem the robust cognitive beamforming design problem. 

In non-CR settings, the study of multiple antenna systems with partial CSI has received considerable 
attention in the past [9], [10]. Specifically, the paper [10] considers the case in which the receiver has 
perfect CSI but the transmitter has only partial CSI (mean feedback or covariance feedback). It was 
proved in [10] that the optimal transmission directions are the same as those of the eigenvectors of 
the channel covariance matrix. However, the optimal power allocation solution was not given in an 
analytical form. A universal optimality condition for beamforming was explored in [11], and quantized 
feedback was studied in [12]. 

In CR settings, power allocation strategies have been developed for multiple access channels (MAC) 
[13] and for point-to-point multiple-input multiple-output (MIMO) channels [14]. Particularly, the 
solution developed in [14] can be viewed as cognitive beamforming since the SU-Tx forms its main 
beam direction with awareness of its interference to the PU. A closed-form method has been present 
in [14]. A water-filling based algorithm is proposed in [13] to obtain the suboptimal power allocation 
strategy. However, the papers [13] and [14] assume that perfect CSI about the link from the SU-Tx to 
the PU is available at the SU-Tx. Due to loose cooperation between the SU and the PU, it could be 
difficult or even infeasible for the SU-Tx to acquire accurate CSI between the SU-Tx to the PU. 

In this paper, we formulate the robust cognitive beamforming design problem as a semi-infinite 
programming (SIP) problem, which is difficult to solve directly. The contribution of this paper can be 
summarized as follows. 

1) Several important properties of the optimal solution of the SIP problem, the rank-1 property, and 
the sufficient and necessary conditions of the optimal solution, are presented. These properties 
would transform the SIP problem into a finite constraint optimization problem. 
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2) Based on these properties, we show that the SIP problem can be transformed into a second order 
cone programming (SOCP) problem, which can be solved via a standard interior point algorithm. 

3) By exploiting the geometric properties of the optimal solution, a closed-form solution for the SIP 
problem is also provided. 

The rest of this paper is organized as follows. Section II describes the SU MISO communication 
system model, and the problem formulation of the robust cognitive beamforming design. Section III 
presents several important lemmas that are used to develop the algorithms. Two different algorithms, 
the SOCP based solution and the analytical solution, are developed in Section V and Section IV, 
respectively. Section VI presents simulation examples, and finally. Section VII concludes the paper. 

The following notation is used in this paper. Boldface upper and lower case letters are used to 
denote matrices and vectors, respectively, (■)^ and (■)^ denote the conjugate transpose and transpose, 
respectively, / denotes an identity matrix, tr(-) denotes the trace operation, and Rank (A) denotes the 
rank of the matrix A. 

II. Signal Model and Problem Formulation 

With reference to Fig. 1, we consider a point-to-point SU MISO communication system, where 
the SU has N transmit antennas and a single receive antenna. The signal model of the SU can be 
represented as y = h^x + n, where y and x axe the received and transmitted signals respectively, 
hg denotes the x 1 channel response from the SU-Tx to the SU-Rx, and n is independent and 
identically distributed (i.i.d.) Gaussian noise with zero mean and unit variance^. Suppose that the PU 
has one receive antenna. The channel response from the SU-Tx to the PU is denoted by an x 1 
vector h. Further, assume that the SU-Tx has perfect CSI for its own link, i.e., hs is perfectly known at 
the SU-Tx. However, due to the loose cooperation between the SU and the PU, only partial CSI about 
h is assumed to be available at the SU-Tx. We assume that ho and R are the mean and covariance 

'since the SU receiver cannot differentiate the interference from the PU from the background noise, the term n can be viewed as the 
summation of the interference and the noise. The variance of n does not influence the algoritfmis discussed later. Moreover, the variance 
of n can be measured at the SU receiver [13]. 
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of h, respectively^. In previous work [10], [15]-[17], partial CSI has been considered in two extreme 
cases in a non-CR setting. One is the mean feedback case, R = a^/, where cr^ can be viewed as the 
variance of the estimation error; and the other is the covariance feedback case, where is a zero 
vector. In this paper, we study the case where the SU-Tx knows both the mean and covariance of h 
in a CR setting. 

The objective of this paper is to determine the optimal transmit signal covariance matrix such that 
the information rate of the SU link is maximized while the QoS of the PU is guaranteed under a robust 
design scenario, i.e., the instantaneous interference power for the PU should remain below a given 
threshold for all the h in the uncertain region. Mathematically, the problem is formulated as follows: 

Robust design problem (PI) : max log(l + hfSh,) 

S>o 

subject to : tr(5) < P, and h^Sh < Pt for {h - ho)^R-\h - ho) < e, 
where S is the transmit signal covariance matrix, P is the transmit power budget, Pt is the interference 
threshold of the PU, and e is a positive constant. The parameter e characterizes the uncertainty of h 
at the SU. According to the definition of the uncertainty in [18], PI belongs to a type of ellipsoid 
uncertainty problem, i.e., the uncertain parameter h is confined in a range of an ellipsoid H, where 
H : {h\{h — Hq)^ R~^{h — Hq) < e}. Thus, the optimal solution of problem PI can guarantee the 
interference power constraint of the PU for all the h E H, and thus the robustness of PI is in the 
worst case sense [19], i.e., in the worst case channel realization, the interference constraint should also 
be satisfied. If the primary transmission does not exist, then the interference constraint is excluded, and 
thus the problem reduces to a trivial beamforming problem. Hence, we only focus on the case where 
the both PU and SU transmission exist. 

Remark 1: An important observation is that the objective function in problem PI remains invariant 
when hg undergoes an arbitrary phase rotation. Without loss of generality, we assume, in the sequel, 
that hs and ho have the same phase, i.e., lm{h^ho} = 0. 

^Due to the cognitive property, we assume that the SU can obtain the pilot signal from the PU, and thus can detect the channel 
information from the PU to the SU. Moreover, since the SU shares the same spectrum with the PU, based on the channel from the PU 
to the SU, the statistics of the channel from the SU to the PU can be obtained [15]. Therefore, we can assume that ho and R are known 
to the SU. 
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Since problem PI has a finite number of decision variable S, and is subjected to an infinite number 
of constraints with respect to the compact set H, problem PI is an SIP problem [20]. One obvious 
approach for an SIP problem is to transform it into a finite constraint problem. However, there is no 
universal algorithm to determine the equivalent finite constraints such that the transformed problem has 
the same solution as the original SIP problem. In the following section, we first study several important 
properties of problem PI, which would be used to transform the SIP problem into its equivalent finite 
constraint counterpart. 

III. Properties of The Optimal Solution 

The maximization problem PI is a convex optimization problem, and thus has a unique optimal 
solution. The following lemma presents a key property of the optimal solution of problem PI (see 
Appendix A for the proof). 

Lemma 1: The optimal covariance matrix S for problem PI is a rank-1 matrix. 

Remark 2: Lemma 1 indicates that beamforming is the optimal transmission strategy for problem 
PI, and the optimal transmit covariance matrix can be expressed as Sopt = Popt'i'opti'^t' where p^pt is 
the optimal transmit power and Vopt is the optimal beamforming vector with ||t'opt|| = 1- Therefore, the 
ultimate objective of problem PI is to determine popt and Vopt- 

According to Lemma 1, a necessary and sufficient condition for the optimal solution of problem PI 
is presented as follows (refer to Appendix B for the proof). 

Lemma 2: A necessary and sufficient condition for 5opt to be the globally optimal solution of problem 
PI is that there exists an hgpt such that 

'-'opt = argmax log(l + hfShs), subject to : tr(5) < p, < p < P, h^^^Sh^p, < Pt, (2) 

S,p 

where 

hopt = argmaxh^ Sopth, for {h — h^)^ R~^{h — Hq) < e. (3) 
h 

Remark 3: The vector h,opt is a key element for all h : (/i — ho)^-R^^(^ — ^o) ^ in the sense that, 
for the optimal solution, the constraint h^^Shopt < Pt dominates the whole interference constraints. 
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i.e., all the other interference constraints are inactive. Thus, if we can determine h^pt, the SIP problem 
PI is transformed into a finite constraint problem (2). It is worth noting that the problem (2) has the 
same form as the problem discuss in [14], in which the CSI on the link of the SU and the link between 
SU-Tx and PU are perfectly known at the SU-Tx. However, unlike the problem in [14], /lopt in (2) is 
an unknown parameter. 

In the following lemma (see Appendix C for the proof), the optimal beamforming vector ?;opt is 
shown to lie in a two-dimensional (2-D) space spanned by Hq and the projection of hg into the null 
space of ho. Define h = h,o/||^o|| and ifi± = where h± = hs — {h hs)h. Hence, we have 

hs = ah,h + bh,h± with ah,,bh, E R. 

Lemma 3: The optimal beamforming vector v^pt is of the form Uyh + byh± with a„, E M. 

Remark 4: According to Lemma 3, we can search for the optimal beamforming vector Vopt on the 
2-D space spanned by h and h±, which simplifies the search process significantly. The optimal v^pt 
found in this 2-D space, is also the globally optimal solution of the original problem PI. As depicted 
in Fig. 2, problem PI is transformed into the problem of determining the beamforming vector Vopt 
in the 2-D space and the corresponding power popt- Combining Lemma 2 and Lemma 3, it is easy to 
conclude that ^opt lies in the space spanned by h and h±. 

IV. Second Order Cone Programming Solution 

In this section, we solve problem PI via a standard interior point algorithm [19], [21], [22]. We 
first transform the SIP problem into a finite constraint problem, and further transform it into a standard 
SOCP form, which can be solved by using a standard software package such as SeDuMi [23]. One 
key observation is that if max^g^^^^ Sh < Pf, i.e., the worst case interference constraint of PI 
is satisfied, then the interference constraint of PI holds. Combining this observation with Lemma 1, 
problem PI can be transformed as: 

Equivalent problem (P2): max 

p>0,\\v\\=l 

(4) 

subject to : p < P, max ph^vv^h < Pt, 
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where 7i(e) := {h\h = Hq + hi}. It is clear that maximizing log(l + ph^vv^hg) is equivalent to 
maximizing \y/ph^v\. By defining w = ^/pv, the objective function can be rewritten as \hfw\. 
Similarly, the interference power can be expressed as \h^w\'^. Thus, problem P2 can be further 
transformed to 

max \h^w\ 
w 

„ „ /- , H , r- 

subject to : \\w\\ < VP, max 

According to the definition of 'H{e), we can rewrite the worst-case constraint in (5) as 

max \h^w\ = max \{hQ + hi)^w\ < \f^f, (6) 

where 7ii(e) := {hi\hi R^^hi < e}. By applying the triangle inequality and the fact that -ye||Qiu|| = 
max\hfw\ for hi E ?ii(e) (refer to Appendix D for details), the interference power can be transformed 
as follows: 

\{ho + hi)^w\ <\h^w\ + \h^w\<\h^w\ + y/^\\Qw\\, (7) 

where Q = A"^/^L/ with A and U being obtained by the eigenvalue decomposition of as 
= U"AU. Moreover, since the arbitrary phase rotation of w does not change the value of the 
objective function or the constraints, according to Remark 1 and Lemma 3, we can assume that w, hg, 
and ho have the same phase, i.e., 

Re{w^hs} > 0, lm{w^ho} = 0, and lm{w^hs} = 0. (8) 

Hence, the interference constraint can be transformed into two second order cone inequalities as follows 

Ve\\Qw\\ + h^w < >/Pt, and ^/^\\Qw\\ - h^w < >/Pf (9) 
By combining (5), (9), with (8), problem PI is transformed into the standard SOCP problem as follows 

max h^w 

(10) 

subject to : \\w\\ < \/p, lm{w^ho} = 0, ^/e\\Qw\\ + hQW < a/^, ^/e\\Qw\\ - h^w < \fPf 
Since the parameters hs and /iq, and the variable w in (10) have complex values, we first convert them to 
its corresponding real-valued form in order to simplify the solution. Define w := [Re{w}^ , lm{w}^]^. 
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h_o := [Re{hoV,lm{ho}T, := [Ro{h,}^ ,lm{h,}T , k := [Im{^o}^, -Re{^}^]^, and Q 
Re{Q} -Im{Q} 

Im{Q} Re{Q} 

We then can rewrife the standard SOCP problem (10) as 



max w 
w 

subject to : WwW < y P, h^w = y/e\\Qw\\ + Hq w < \/ Pt, \^\\Qw\\ — w < \/ Pt- 



(11) 



Problem (11) can be solved by a standard interior point program SeDuMi [23], which has a polyno- 
mial complexity. In the next section, we develop an analytical algorithm to solve problem PI, which 
reduces the complexity of the interior point based algorithm substantially. 

V. An Analytical Solution 

In this section, we present a geometric approach to problem PI. We begin by studying a special case, 
the mean feedback case, i.e., -R = cr^/. Due to its special geometric structure, the mean feedback case 
problem can be solved via a closed-form algorithm. We next show that problem PI can be transformed 
into an optimization problem similar to the mean feedback case. Based on the closed-form solution 
derived for the mean feedback case, the analytical solution to problem PI with a general form of a 
covariance matrix R is presented in Subsection V-B. 

A. Mean Feedback Case 

Based on the observation in Lemma 1 and the definition of the mean feedback, the special case of 
problem PI with mean feedback can be written as follows. 

Mean feedback problem (P3): max log(l +phfvv^hs) 

p>om=i (12) 

subject to : p < P, ph^vv^h < Pt, for \\h - ho\\^ < ecr^. 
Problem P3 has two constraints, i.e., the transmit power constraint and the interference constraint. 
Similar to the idea in [13], the two-constraint problem is decoupled into two single-constraint subprob- 
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lems: 



Subproblem 1 (SPl): 



max 

P>0,||ti|l=l 

subject to : 



log(l + phg vv 




(13) 



p< P. 



(14) 



Subproblem 2 (SP2): 



p>o,|ii;|j=i 
subject to : 



max 



log(l +ph^ vv 




(15) 



ph 



vv"h < Pt, for \\h-hof < ea\ 



(16) 



In the sequel, we present the algorithm to obtain the optimal power popt and the optimal beamforming 
vector Vopt for both subproblems in subsection V-A.l, and describe the relationship between the 
subproblems and problem P3 in subsection V-A.2. 

1) Solution to subproblems: For SPl, the optimal power is constrained by the transmit power 
constraint, and thus popt = P- Moreover, since there does not exist any constraints on the beamforming 
direction, it is obvious that the optimal beamforming direction is equal to hs, i.e., Vopt = ^s/||^s||. 
Thus, the optimal covariance matrix Sopt for SPl is Phsh^ /WhsW^. In the following, we focus on the 
solution to SP2. 

SP2 has infinitely many interference constraints, and thus is an SIP problem too. By following a 
similar line of thinking as in Lemma 2, SP2 can be transformed into an equivalent problem that has 
finite constraints (refer to Appendix E for the proof) as follows. 

Lemma 4: SP2 and the following optimization problem: 



where /lopt = + V^crv, have the same optimal solution. 

According to Lemma 4, problem (17) has the same optimal solution as SP2. Moreover, according 
to Lemma 3, the optimal solution v of problem (17) lies in the plane spanned by h and h±. We next 
apply a geometric approach to search the optimal solution, i.e., by restricting our search space to a 2-D 
space. As shown in Fig. 3, assume that the angle between v and ho is /3, and the angle between hg 
and ho is a. It is easy to observe that < a < Since v lies in a 2-D space, v can be uniquely 

^This follows because if a > tt/2, we can always replace hs by —hs without affecting the final result, and the angle between —hs 
and ho is less than 7r/2. 



max log(l + ph vv 

p>0,||t!||=l 



hs), subject to : ph ^vv h^^^ < Pt 



(17) 
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identified by the angle [3. Hence, we need only to search for the optimal angle /3opt. By exploiting the 
relationship between p, v, and P, the two-variable optimization problem (17) can be further transformed 
into an optimization problem with a single variable P, which can be readily solved. 

By observing Fig. 3, the angle between hg and v is /3 — a, and hence the objective function of (17) 
can be expressed as 

max log(l + phfvv^hs) = maxlog ( 1 + cos^(/3 — a)). (18) 

||'y||=i 13 \ / 

Clearly, the maximum rate is achieved if the following function 

/(/?) := p\\hs\\\os\(3 - a) (19) 

is maximized. 

Moreover, it can be proved by contradiction that the interference constraint is satisfied with equality, 
i.e., h^p^Shopt = Pt- Thus, we have 

ph^^^vv^'K^, = p{ho + V~e(rvfvv''{ho + V^av) = p{\\ho\\ cos/5 + v^a)' = Pf (20) 

Hence, the interference constraint is transformed into 

P^ 

I /loll cos/3 + i/ecr) 

By substituting (21) into (19), we have 



P = ^^. Z (21) 



rfn^ || . ||2 2fn N | /^s P^f COS^ /^ - « 

f{P)=p\\hs\\ COS {f3-a) = —— — —2. (22) 

(||/io||cos(/5) + 

Thus, the optimal /5opt can be expressed as 

Popt = arg max / [p) = arg max . (23) 

{\\ho\\cos{l3) + ^/ea) 

The problem of (23) is a single variable optimization problem. It is easy to observe that the feasible 
region for /? is [a,n/2]. According to the sufficient and necessary condition for the optimal solution of 
an optimization problem, P^pt lies either on the border of the region (a or 7r/2) or on the point which 
satisfies df{p)/dp = 0. Since 

dfl^^-^ 2\\hs\\'^ Pt cos{P - a)^sina - sin(/3 - a)v/eo-/||/io||) 

~W ~ ||/ior(cos/3+v^a/||/io||)' ' ^^"^^ 

January 4, 2009 DRAFT 



IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS (REVISED) 11 

we can obtain a locally optimal solution (3i = sin~^ ^ ||/i-o|[&inQ j _|_^ solving the equation df{(3)/d(3 = 
0. In the case when ll^ojL^"^" > l, is a non-decreasing function. Hence, the optimal [3 is 7r/2, and 
we define = —oo for this case. Therefore, the globally optimal solution is 

/3opt = argmax(/(a),/(7r/2),/(A)). (25) 

The optimal power popt can be further obtained by substituting /3opt into (21). According to the 
definition of [3 and Lemma 3, we have 

i^opt = ttyh + b^h±, (26) 

where = cos(/3opt) and by = sin(/3opt)- In summary, SP2 can be solved by Algorithm 1 as described 
in Table I. 

2) Optimal solution to problem P3.' In the preceding subsection, we presented the optimal solutions 
for the two subproblems. We now turn our attention to the relationship between problem P3 and the 
subproblems, and present the complete algorithm to solve problem P3. Since the convex optimization 
problem P3 has two constraints, the optimal solution can be classified into three cases depending 
on the activeness of the constraints: 1) only the transmit power constraint is active; 2) only the 
interference constraint is active; and 3) both constraints are active. Relying on this classification, the 
relationship between the solutions of problem P3 and the two subproblems is described as follows 
(refer to Appendix F for the proof). 

Lemma 5: If the optimal solution Si of SPl satisfies the constraint of SP2, then Si is the optimal 
solution of problem P3. If the optimal solution ^2 of SP2 satisfies the constraint of SPl, then ^2 
is the optimal solution of problem P3. Otherwise, the optimal solution of problem P3 simultaneously 
satisfies the transmit power constraint and h^^^Sho^x < Pt with equality. 

Remark 5: To apply Lemma 5, we need to test whether Si and ^2 satisfy both constraints. The 
condition that Si satisfies the interference constraint is 

Pint < Pt, where Pi„t = max h^Sih, for - /iqIP < ea^, (27) 

h 

where Pint can be obtained by the method discussed in Appendix D. The condition that S2 satisfies 
the transmit power constraint is tr(52) < P. 
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We next discuss the method for finding the solution in the case where neither Si nor S2 is the 
optimal solution of problem P3. Similarly to the method in the preceding subsection, we solve this 
case from a geometric perspective. According to Lemma 5, in the case in which neither Si nor ^2 is 
the feasible solution, the optimal covariance 5opt must satisfy both constraints with equality, i.e., 

Popt = P, and Popt^5t'"°pt'"St'*°pt = ^t- (28) 
Combining these two equalities, we have 

P{\\ho\\cos{(3) + ^/~eay = Pt. (29) 

Thus, 

/3„p, = arccos(^^^^/5^). (30) 

Based on P^pu we can obtain v^pt from (26). We summarize the procedure called Algorithm 2, which 
solves the case where both constraints are active for problem P3, in Table II. Furthermore, we are 
now ready to present the complete algorithm, namely Algorithm 3, to solve problem P3 in Table III. 

In Algorithm 3, we obtain the optimal solutions to SPl and SP2 and the optimal solution to the 
case where both constraints are active separately. According to Lemma 5, the final solution obtained 
in Algorithm 3 is thus the optimal solution of problem P3. 

Proposition 1: Algorithm 3 obtains the optimal solution of problem P3. 

B. The Analytical Method for Problem P 1 

In the preceding subsection, the mean feedback problem P3 is solved via a closed-form algorithm. 
Unlike problem P3, problem PI has a non-identity-matrix covariance feedback. To exploit the closed- 
form algorithm, we first transform problem PI into a problem with the mean feedback form as follows. 

Equivalent problem (P4): maxlog(l +phg vv hg) 

subject to : p|| A^/^i;f < P, ph^vv^h < Pt, for \\h - hoW"^ < e, 
where R^^ := AU obtained by eigen-decomposing R^^, h := A^^'^Uh, ho := A^/'^UHq, 
hg := A^/'^Uhs, and v := A~^/'^Uv. By substituting these definitions into (31), it can be observed that 
the achieved rates and constraints of both problem PI and P4 are equivalent. Thus, the optimal solution 
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of PI can be obtained by solving its equivalent problem P4. Moreover, the optimal beamforming vector 
■Uopt of problem P4 can be easily transformed into the optimal solution Vopt for problem PI by letting 
Vopt = A^^'^Vopf Note that it is not necessary that = 1 in (31). 

In the preceding subsection, decoupling the multiple constraint problem into several single constraint 
subproblems facilitates the analysis and simplifies the process of solving the problem. For problem P4, 
it can also be decoupled into two subproblems as follows. 

Subproblem 3 (SP3): maxlog(l + phf vv" h,) (32) 

subject to : p||A^/2^||^ < p. (33) 
Subproblem 4 (SP4): max\og{l + phf vv^ hs) (34) 

p,V 

subject to : ph^vv^h < Pt for \\h - hof < e. (35) 

It is easy to observe that SP3 is equivalent to SPl, and the optimal transmit covariance matrix of 
SP3 can be obtained in the same way as that for SPl. Moreover, SP4 is the same as SP2, and thus 
it can be solved by Algorithm 1 discussed in Subsection V-A.l. 

The relationship between problem P4 and subproblems SP3 and SP4 is similar to the one between 
P3 and corresponding subproblems as depicted in Lemma 5, i.e., if either optimal solution of SP3 or 
SP4 satisfies both constraints, then it is the globally optimal solution; otherwise, the optimal solution 
satisfies both constraints with equalities. We hereafter need to consider only the case in which the 
solutions of both subproblems are not feasible for problem P4. For this case, the two equality constraints 
can be written as follows. 

||A^/^t;|| = 1, and max (h^vv^h) = for ||^ - hof < e. (36) 

Assume that the angle between ho and v is (3, and that p := II^H. Similar to Lemma 3, the optimal v 
lies in a plane spanned by h and h±, where h = hQ/\\ho\\, h± = and h± = hg — {h^hs)h. 

Thus, if we can determine /3 and p from (36), then the optimal v can be identified by 

V = p{cos{(3)h + sm.{l3)hA_). (37) 
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Based on the new variables l3 and p, the constraints (36) can be transformed as follows. 



P 



cos(/5)/i + sm{(3)hj_) = 1, 



(38) 



and, p(cos(^)||/io|| + V^) 




(39) 



According to (38), we have 



1 



P = 



(40) 



AV2(cos(/?)^ + sin(/3)/ix) 



Substituting (40) into (39), we have 




cos(^)||^o|| + v^- 



(41) 



Hence, the optimal f3 can be obtained by solving (41), and v^pt can be obtained by substituting /3 
into (37). In summary, the procedure to solve the case in which both constraints are active is listed 
as Algorithm 4 in Table IV. Moreover, we are now ready to present the complete algorithm, namely 
Algorithm 5, for solving problem PI in Table V. 

In Algorithm 5, we obtain the optimal solutions to SP3 and SP4 and the optimal solution to the 
case where both constraints are active separately. According to Lemma 5, the final result obtained in 
Algorithm 5 is thus the optimal solution of problem PI. 

Proposition 2: Algorithm 5 achieves the optimal solution of problem PI. 

Remark 6: The complexity of the interior point algorithm for the SOCP problem (1 1) is 0{N^'^ log(i)), 
where e denotes the error tolerance. For Algorithm 5, a maximum of C(log(i)) operations is needed to 
solve (41), and the complexity for each operation is C(log(A^^)). Hence, the computation complexity 
required for Algorithm 5 is 0{N'^ \og(-)), which is much less than that of the interior point algorithm. 



Computer simulations are provided in this section to evaluate the performance of the proposed 
algorithms. In the simulations, it is assumed that the entries of the channel vectors hg and ho are 
modeled as independent circularly symmetric complex Gaussian random variables with zero mean and 
unit variance. Moreover, we denote by /i the distance between the SU-Tx and the SU-Rx, and by I2 



VI. Simulations 
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the distance between the SU-Tx and the PU. It is assumed that the same path loss model is used to 
describe the transmissions from the SU-Tx to the SU-Rx and to the PU, and the path loss exponent is 
chosen to be 4. The noise power is chosen to be 1, and the transmit power and interference power are 
defined in dB relative to the noise power. For all cases, we choose Pt = dB. 

A. Comparison of the Analytical Solution and the Solution Obtained by the SOCP Algorithm 

In this simulation, we compare the two results obtained by a standard SOCP algorithm (SeDuMi) 
and Algorithm 3. We consider the system with = 3, /2//1 = 2, and P ranging from 3 dB to 10 dB. 
In Fig. 4, we can see that the results obtained by different algorithms coincide. This is because both 
algorithms determine the optimal solution. Compared with the SOCP algorithm solution. Algorithm 3 
obtains the solution directly, and thus it has lower complexity. In Fig. 5, we compare the two results 
obtained by SeDuMi and Algorithm 5. We consider the system with A^ = 3, P = 5 dB, and I2/I1 
ranging from 1 to 10. The covariance matrix R is generated by Ri, where each element of Ri 
follows Gaussian distribution with zero mean and unit variance. From Fig. 5, we can see that the results 
obtained by the two algorithms coincide again. Moreover, we note that the achievable rate with e = 0.2 
is always greater than or equal to the rate with e = 0.3, since a larger e corresponds to the stricter 
constraints. 

B. Ejfectiveness of the Interference Constraint 

In this simulation, we apply Algorithm 3 to solve problem P3. In Fig. 6, we depict the achievable 
rate versus the ratio h/h under different transmit power constraints. The increase of the ratio h/h 
corresponds the decrease of the interference power constraint. As shown in Fig. 6, with an increase of 
I2/I1, the achievable rate increases due to the lower interference constraint. Until the ratio I2/I1 reaches 
a certain value, the achievable rate remains unchanged, since the transmit power constraint dominates 
the result, and the interference constraint becomes inactive. 

C. The Activeness of the Constraints 

In this simulation, we compare the achieved rates of problem PI with a single transmit power 
constraint, a single interference constraint and both constraints. Here, we choose A^ = 3, e = 0.2, and 
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generate R in the same way as in the first simulation example. Fig. 7 plots three achievable rates for 
different constraints, respectively. It can be observed from Fig. 7 that the rate under two constraints is 
always less than or equal to the rate under a single constraint. Obviously, this is due to the fact that 
extra constraints reduce the degree of freedom of the transmitter. 

VII. Conclusions 

In this paper, the robust cognitive beamforming design problem has been investigated, for the SU 
MISO communication system in which only partial CSI of the link from the SU-Tx to the PU is available 
at the SU-Tx. The problem can be formulated as an SIP optimization problem. Two approaches have 
been proposed to obtain the optimal solution of the problem; one approach is based on a standard 
interior point algorithm, while the other approach solves the problem analytically. Simulation examples 
have been used to present a comparison of the two approaches as well as to study the effectiveness 
and activeness of imposed constraints. 

This work initiates research in robust design of cognitive radios. We are currently extending these 
methods to the more general case with multiple receive antennas and multiple PUs. Other interesting 
extensions include more practical scenarios, such as the case in which the SU channel information is 
also partially known at the SU-Tx. 

Appendix 

A. Proof of Lemma 1: Problem PI involves infinitely many constraints. Denote the set of active 
constraints by C, the cardinality of the set C by K, and the channel response related to the kth element 
of the set C by hk. According to the Karush-Kuhn-Tucker (KKT) conditions for PI, we have: 

K 

hsil + h^Shs)-'h^ + $ = A/ + 5^ fiih,h^, (42) 

i=l 

tr(*5) = 0, (43) 

where $ is the dual variable associated with the constraint S > 0, and A and /ij are the dual variables 
associated with the transmit power constraint and the interference constraint, respectively. First, we 
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assume that A 7^ 0, and thus the rank of the right hand side of (42) is N. Since the first term on the 
left hand side of (42) has rank one, we have 

Rank($) > TV - 1. (44) 

Moreover, since 5 > and $ > 0, from (43) we have tr($5) = tr{U^\US) = tr(At/5C/^) = 
tr(AS') = 0, where AU is the eigenvalue decomposition of matrix and S := USU^ . By 
applying eigenvalue decomposition to S, we have S := Yli'^i^i^? ^ where is the ith eigenvalue 
and Si is the corresponding eigenvector. We next show Rank(S') + Rank(#) < by contradiction. 
Suppose that Rank(S') + Rank($) > N. Then, there exists an index j such that the jth element of Sj 
and the jth diagonal element of A are non-zero simultaneously. Thus, it is impossible that the equation 
tr(AS') = holds. It follows that Rank(5) + Rank($) < A^. Combining this with (44), we have 
Rank(5) < 1. 

Second, we assume that A = in (42). In this case, S must lie in the space spanned by hi, 
i = 1, ■ ■ ■ , K. Let the dimensionality of the space be M. Therefore, we can restrict Rank($) < M. 
Thus, the reminder of the proof is the same as that of the case A 7^ 0, and the proof is complete. ■ 

B. Proof of Lemma 2 : First, we consider the sufficiency part of this lemma. We assume that there 
exists a covariance matrix Sopt and an ^opt that satisfy the conditions (2) and (3) simultaneously. 
Since Sopt satisfies both the transmit power constraint and the interference constraint, Sopt is a feasible 
solution for problem PI. Moreover, if we assume that there exists another solution Sg, which results in 
a larger achievable rate for the SU link, then a contradiction will be derived. Without loss of generality, 
we assume that the constraint set, which consists of all the active interference constraints for Ss, is 
denoted by T. We divide the set T into two types: one type is /lopt £ and the other type is hopi ^ T. 

Assume that Cg and Copt are the achievable rates for the covariance matrices Sg and Sopt, respectively. 
In the case of hopt G T, we have Cg < Copt, since Copt is obtained with fewer constraints. Since problem 
PI is a convex optimization problem that has a unique optimal solution, Sopt is indeed the optimal 
solution. In the case of /lopt ^ we can observe that Sopt satisfies the constraints in T, and Sg satisfies 
the constraint /lopt- According to the lemma in [13], this case does not exist. 

We next proceed to prove the necessity part. Suppose that Sopt is the optimal solution of problem 
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PI. According to Lemma 1, we have Sopt = Popt'^optf^f Thus, problem PI is equivalent to 



opt" 

max log(l + hfShs) 

S>o (45^ 



subject to : tr(5) < Popt, h^Sh < Pt, for (h - ho)^R-\h - ho) < e. 
According to Lemma 6, there is a unique 



e 



^opt = ho + W— aRVopt, (46) 

V ''^opt-K'i^opt 

which is the optimal solution of max/^g^^^-, Sh < Pt. Thus, for problem (45), only tr(S') < Popt 
and h^^^Shopf < Pt are active constraints. Thus, it is obvious that problem (45) and problem (2) have 
the same optimal solution. Hence, the proof is complete. ■ 

C. Proof of Lemma 3 : The proof of Lemma 3 is divided into two parts. The first part is to prove 
that Vopi is in the form of ai,h + Pi,h±, where E C and E C. The second part is to prove a„ E M 
and Pv G In the following proof, we assume that E C are some proper complex scalars. 
According to Lemma 2, and Theorem 2 in [14], we have 

Vopt = tti^opt + <y2hs- (47) 

According to Lemma 6, we have 

^opt = ho + asVopt = ho + a3{aihopt + a2hs) = ho + aia3hopt + 0203/15- (48) 

According to (48), it can be observed that /lopt can be expressed by the linear combination of ho and 
hg, where the coefficients are complex. Combining this with (47), we have Vopt = O!iho + a^hg, where 
04 G C and 0:5 E C. Moreover, since both ho and hg can be expressed as a linear combination of h 
and h±, we have i^opt = ayh + Pyh±. Since rotating v^pt does not affect the final result, we can assume 
a„ E M. 

We next prove that /3„ G M by contradiction. At first, we assume that Py = a + jb ^ M. Then we 
can find an equivalent (3y = y/a"^ + b'^ E M which is a better solution of problem PI than Py. Assume 
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that Vopt = C(vh + i3vh±. It is clear that \\Vopt\\ = ll'^'optllj and the interference caused by Vopt is 



P'*St'^°pt*St'*opt =p{ho + / ai?i;opt)^'t'opt*^t(^o + , / Q-R^'opt) (49) 



'^opt-"''^opt y t^opt-"''^opt 



=p(a^\\ho\\ + / T-H"^7T — a^vfptRvoptY, (50) 
which is equal to that of Vopt- However, the corresponding objective function with Vopt is 
log(l + phfvoptV^pAs) = log(l + p{ah,h + bh,h±)" {a^h + (3^h±){ayh + l3yh±)^{ah^h + bh^h±)) 

= log{l +p{ah,a^ + bhjv){ah,a^ + hsPv)), (51) 
and the objective value with Vopt is 

log(l + phfvoptv^pihs) = log(l + p{ah,h + bh,h±)" {a^h + l3^h±){ayh + l3yh±)"{ah,h + bh,h±)) 

= log{l +p{ah,ay + bhjv){ah,(y^ + bhjif)). (52) 
According to (51) and (52), we can conclude that Vopt is a better solution. The proof follows. 

D. Lemma 6 and its proof: 
Lemma 6: For the problem 

max ph^vv^h, subject to: {h — Hq)^ R~^{h — Hq) < e, (53) 
h 

where p, v, and Hq are constant, the optimal solution is 



^max = ho + ^ ^J'^^ aRv, where a = v^ho/\v"ho\. (54) 
Proof: The objective function ph^vv^h is a convex function. The duality gap for a convex 
maximization problem is zero. The Lagrangian function is 

L{h,X) =ph"vv^h- x(^{h-ho)"'R-\h-ho) -e^ (55) 

where A is the Lagrange multiplier. According to the KKT condition, we have ^ = 2pvv^h — 
2XR^\h - ho) = 0. Thus, 

p{v"h)v = XR-\h - ho). (56) 
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We have /imax = + baRv, where 6 G M, a G C, and |a| = 1. Since (h — hQ)^R ^{h — h 



we have h = \^/V v^R^v. Moreover, by observing (56), we have a = tv^h = tv^{ho + baRv) = 
tv^fiQ + tbav^Rv, where t is a real scalar such that \tv^h\ = 1. Thus, we have v^Hq/Iv^Hq] = a. 
The proof follows immediately. ■ 

E. Proof of Lemma 4: Similar to the proof of Lemma 2, we can show that the problem 

5opt = argmaxlog(l + hfSh^) subject to : h^^^Shopt < Pt, (57) 

S,p 

where hopt = argmax^ Sopth, for {h — h^)^ R^^ {h — Hq) < e, is equivalent to SP2. 

Since Sopt is a rank-1 matrix, according to Lemma 6, we have /iopt = + ^/ecrv. Combining this 
with (57), we have 5opt = argmax^plog(l + Shg) s.t. : (/iq + y/eav)^ S{ho + y/eav) < Pt, 
which is equivalent to (17). The proof is complete. ■. 

F. Proof of Lemma 5: Assume that Sopt is the optimal solution for problem P3. If Si satisfies the 
interference constraint, then Si is a feasible solution for problem P3. The optimal rate achieved by 
Sopt cannot be larger than that of Si, since the constraint of SPl is a subset of problem P3. Similarly, 
we can prove the second part of the Lemma. We now focus on the third part of this lemma. For problem 
P3, at least one of tr(S') < P and /i^jShopt < Pt is an active constraint, since if neither of them 
is active, we can always find an e such that Sopt + el is a feasible and better solution. Moreover, if 
only tr(S') < P is active, then Si is the optimal solution, which contradicts with /i^jSi/iopt > Pt- 
Similarly, it is impossible that only h^^^Shopt < Pt is active. Therefore, both constraints are active 
constraints. ■ 
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TABLE I 
The algorithm for SP2. 



Algorithm 1 

1. Compute Popt through (25), 

2. Compute popt according to (21), 

3. Compute Vopt according to (26), 

4. <S'opt = Popt'l^opt^^i^f 



TABLE II 

The ALGORITHM FOR PROBLEM P3 IN THE CASE WHERE TWO CONSTRAINTS ARE SATISFIED SIMULTANEOUSLY. 



Algorithm 2 

1. Compute /3opt through (30), 

2. Based on (26), compute Uopi, 

3. tS'opt ~ -P^opt^opt- 



TABLE III 

The complete algorithm for problem P3. 



Algorithm 3 

1. Compute the optimal solution Si — Phsh^ /\\hs\\^ for SPl, 

2. Compute the optimal solution S2 for SP2 via Algorithm 1, 

3. If Si satisfies the interference constraint, then Si is the optimal solution, 

4. Elsif •S'2 satisfies the transmit power constraint, then S2 is the optimal solution, 

5. Otherwise compute the optimal solution via Algorithm 2. 
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TABLE IV 

The algorithm for problem P4 in the case where two constraints are satisfied simultaneously. 
Algorithm 4 

1. Compute P via (41), and compute v via (37), 

2. Based on the relationship between v and v, compute Vopt, 

3. tS'opt = -P'^'opt'^^opt- 



TABLE V 

The complete algorithm for problem PI. 

Algorithm 5 

1. Compute the optimal solution S3 = Phsh" /\\hs\\^ for SP3, 

2. Compute the optimal solution •S'4 for SP4 via Algorithm 4, 

3. If Ss satisfies the interference constraint, then S3 is the optimal solution, 

4. Elsif •S'4 satisfies the transmit power constraint, then S'4 is the optimal solution, 

5. Otherwise compute the optimal solution through Algorithm 4. 




SU-Tx 



Fig. 1. The system model for the MISO SU network coexisting with one PU. 
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Fig. 2. The geometric explanation of Lemma 3. The ellipse is the projection of h := {{h — ho)^ R ^{h — ho) = e} on the plane 
spanned by h and h±. 



January 4, 2009 



DRAFT 



IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS (REVISED) 



25 




January 4, 2009 



DRAFT 



IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS (REVISED) 



26 




Fig. 4. Comparison of thie results obtained by the SOCP algorithm and Algorithm 3. 



January 4, 2009 



DRAFT 



IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS (REVISED) 



27 




Fig. 5. Comparison of thie results obtained by the SOCP algorithm and Algorithm 5. 
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Fig. 6. Effect of i2/h on the achievable rate of the CR network (e = 1, = 3). 
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Fig. 7. Comparison of tlie rate under different constraints of problem PI. (i) tlie maximal rate subject to interference constraint and 
transmit power constraint simultaneously; (ii) the maximal rate subject to a single transmit power constraint; (iii) the maximal rate subject 
to a single interference constraint. 
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